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Abstract 

Specific heat and muon spin rotation and relaxation data are presented for two 
geometrically frustrated systems: Gd2Ti207 which antiferromagnetically orders and 
Yb2Ti207 which presents dynamical short range correlations at low temperature. 
The muon data help to characterize the spin dynamics of these two compounds. 
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Geometrically derived magnetic frustration arises when the spatial arrange- 
ment of the spins is such that it prevents the simultaneous minimization of all 
interaction energies [1]. Compounds of interest for the investigation of frus- 
tration are, for example, the crystallographically ordered pyrochlore structure 
compounds R2T12O7 where the rare earth ions (R) form a sub-lattice of corner 
sharing tetrahedra. 

We shall focus on two pyrochlore compounds Gd 2 Ti 2 07 and Yb 2 Ti 2 7 . Re- 
ports of experimental data concerning these two systems were already pub- 
lished by some of us, see Refs. [2,3,4,5]. A summary of the muon spin rotation 
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Fig. 1. Thermal dependence of the specific heat divided by the temperature mea- 
sured for a crystal of Gd2Ti207. The data are presented per mole of Gd. The 
maxima of the specific heat at T c \ = 1.02 K and T C 2 = 0.74 K are the signatures of 
the magnetic phase transitions. The insert displays the computed entropy variation 
as a function of the temperature between ~ 85 mK and 5 K. 



and relaxation (/xSR) data has been recently provided [6]. Here we report un- 
published data obtained from specific heat and /zSR measurements. The latter 
measurements were either performed at the ISIS /iSR facility in UK or at the 
Swiss Muon Source of the Paul Scherrer Institute in Switzerland. We stress 
that the specific heat measurements are essential to test the sample quality 
in geometrically frustrated systems. For example, in the case of Gd 2 Ti 2 07, 
the heat treatment affects the specific heat data in two ways. First it notably 
increases the intensity of the signal and second it slightly decreases the tem- 
perature of the two peaks corresponding to the magnetic transitions. As far 
as fiSR spectra are concerned, the values for the measured local fields are 
noticeably changed by the heat treatment. 

In a first step we consider Gd2Ti 2 0y. Figure 1 shows the specific heat mea- 
surements and the variation of entropy which is deduced from these data. The 
entropy variation up to ~ 5 K is only 90 % of its expected value for a spin 7/2 
system. This missing entropy is probably related to the persistent spin dy- 
namics observed in this system deep in the ordered state [2]. The specific heat 
structure observed at T c \ ~ 1.02 K for Gd2Ti 2 07 corresponds to a long-range 
magnetic ordering of the Gd magnetic moments as shown by the spontaneous 
muon spin precession detected below T cl in zero-field. In fact, we observe two 
Bessel-like oscillations, each characterized by the maximum field at the muon 
site, -B max ; see Fig. 2. Note that the transition at T c2 ~ 0.74 K is hardly vis- 
ible in the temperature dependence of B max . This is not surprising since the 
difference between the magnetic structures adopted by Gd2Ti 2 0y above and 
below T C 2 is small: the change concerns only a quarter of the Gd 3+ moments 
which are disordered above T c2 and very weakly ordered below [7]. 
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Fig. 2. Temperature dependence of -B max measured for Gd2Ti20y (two fields per 
temperature point can be resolved up to 0.6 K and only one field at higher temper- 
atures) and the associated exponential damping rate Ax for the higher B mSLX value. 
The damping rate is negligible for the smaller -B max - The line is the result of a fit for 
the larger field with B max (T) oc [1 - (T/T cl ) a f, setting a = 2 and (3 = 0.365 as ex- 
pected for a 3D Heisenberg antiferromagnet characterized by a second order phase 
transition at T = T c \. The fit is numerically reasonable and physically justified since 
-Bmax is expected to track the order parameter. The contact field term, which might 
have a temperature dependence different from that of the order parameter, should 
indeed be negligible for an insulating compound such as Gd2Ti207. 

Although the thermal behavior of the spin-lattice relaxation, Xz, at high tem- 
perature is not directly related to the frustrated nature of the compounds of 
interest, it is still of interest. We present A^(T) above T cl for Gd 2 Ti 2 7 in 
Fig. 3. When the temperature is decreased from room temperature towards 
T c i, ^z{T) first slightly decreases and reaches a local minimum around 10 K, 
before the slowing down of the spin dynamics causes a marked increase of 
\z{T) when approaching T cl . The initial slight decrease is the signature of the 
antiferromagnetic spin correlations which build up in Gd 2 Ti 2 7 as the system 
is cooled down. It is fitted to a model in which the spin correlation function 
is expanded versus (/cbT) -1 [6,3] (see line in Fig. 3). 

A remarkable output of most /iSR measurements for geometrically frustrated 
magnetic materials has been the observation of persistent spin dynamics at 
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Fig. 3. The spin-lattice relaxation rate Xz versus temperature measured in the 
paramagnetic phase for crystals of Gc^T^CV The data have been recorded for 
two orientations of the initial muon polarization, parallel to either a [111] or [110] 
crystallographic axis: they show no difference. The line is the result of a fit with 
\ Z (T) = \ ( z s \l + T /T) where X^ = 1.45 (37) /xs" 1 and T = -2.3(6) K. The 
further increase of Xz as the sample is cooled down below ~ 10 K reflects the slowing 
down of the magnetic fluctuations when approaching T c \ from above. This effect is 
not accounted for by the model. 



low temperature. That dynamics has recently been shown to be the signature 
of a large density of states at low energy characterized by a gap linear in tem- 
perature [2] . The persistent dynamics is usually recognized by the observation 
of a finite value for Xz- This is the case for Gd2Ti207 [2]. If the spin dynamics 
is particularly slow, Pzif) is no longer stretched-like but Kubo-Toyabe-like. 
This type of relaxation function was found in Yb2Ti 2 07 [4] below the tem- 
perature of the specific heat peak shown in Fig. 4. The spin dynamics is then 
revealed by the finite slope of Pz(t) at large time. In order to ascertain the 
interpretation of the Yb2Ti 2 07 spectra, we have investigated their field de- 
pendence. As shown in Fig. 5, the relaxation is quenched by a longitudinal 
field. The fits presented as full lines in this figure are done using the function 
a Pz(t) = gii-PgbgW + a 2 where -Pgbg(^) is the Gaussian-broadened Gaussian 
function [8] and the second term (a 2 ) represents the background signal which 
mainly arises from muons stopped in the silver sample holder. A conventional 
Kubo-Toyabe function instead of -Pgbg(^) does not properly account for the 
data. This is illustrated by the dashed line in Fig. 5 which is the best fit for the 
spectrum recorded in a 2 mT longitudinal using the Kubo-Toyabe function: 
the dashed line curve differs notably from the data for t < 1 /zs. It means that 
the muon spins detect a field distribution influenced by disorder. This proba- 
bly results from the geometrical frustration of the magnetic interactions. We 
also note that a dynamical spin-glass model [9] does not either provide a fit 
of the early times data as good as the -Pgbg(^) function. 

The model behind this Gaussian-broadened Gaussian function assumes that 
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Fig. 4. Specific heat variation for the powder sample of Yb2Ti207 which was used 
in a previously published study [4]. The data are plotted per mole of Yb. The 
maximum of the peak is around 0.25 K, which is the temperature of a sharp change 
in the spin dynamics [4]. The insert shows the entropy variation deduced from the 
specific heat data. 



the muons see different environments, each of them being characterized by 
a dynamical Kubo-Toyabe function associated with a field width A/7^. To 
account for the slight differences between the environments, A is then assumed 
to be Gaussian distributed with a mean- value Ao and a width w. The common 
parameters for the five spectra presented in Fig. 5 are the initial asymmetries 
of the /iSR signals related to the sample and the background, respectively ai 
= 0.170 and a 2 = 0.065, A = 4.9/xs -1 (A /7 M = 5.7 mT), the correlation 
frequency of the field at the muon v = 0.85 MHz and w/A Q = 0.38. Only 
the field values are changed. As expected from the important effect of modest 
fields in the shape of the /xSR spectra, v is found in the megahertz range. 
The remarkable result is that the fits work only if we take the values for the 
field smaller than the actual applied field. We have B meas = 2.5, 5.0, 9.2, 17.2 
and 36 mT instead of B ext = 10, 20, 50, 100 and 200 mT respectively. This 
remarkably feature has already been encountered for a Kagome-like material 
[10]. It has been attributed to the intermittent nature of the spin dynamics. 
In our case the ratio B meas /B ext depends on B ext which would imply that the 
applied field affects the system. We note that Aq/7^ is reduced compared to 
the high temperature field width of 80 mT [11] measured above the specific 
heat anomaly at ~ 0.25 K. This reduction is consistent with the sporadic 
nature of the field at the muon. More work is yet needed for an interpretation 
of the magnitude of this reduction. 
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Fig. 5. Effect of a longitudinal magnetic field on the longitudinal polarization func- 
tion Pz(t) recorded in Yb2Ti207 at low temperature. The spectra are presented in 
absolute scale, i.e. they are not shifted vertically. The different lines are results of 
fits described in the main text. 
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